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LAGRANGIAN ASPECTS OF THE AXISYMMETRIC EULER 

EQUATION 


STEPHEN C. PRESTON AND ALEJANDRO SARRIA 

Abstract. In this paper we are interested in geometric aspects of blowup in the axisym- 
nretric 3D Euler equations with swirl on a cylinder. Writing the equations in Lagrangian 
form for the flow derivative along either the axis or the boundary and imposing oddness on 
the vertical component of the flow, we extend some blowup criteria due to Chae, Constantin, 
and Wu related to assumptions on the sign of the pressure Hessian. In addition we give a 
geometric interpretation of the results, both in terms of the local geometry along trajectories 
and in terms of the Riemannian geometry of the volume-preserving diffeomorphism group. 


1. Introduction 

The question of whether smooth solutions of the three-dimensional Euler equations can 
break down in finite time is a long-standing open problem; see Constantin jBJ for a history 
and survey of results. Even in the axisymmetric case, where the velocity components do 
not depend on the angular coordinate 6 , the question is still open, although if in addition 
the angular velocity is assumed to be zero (that is, axisymmetric flow without swirl), global 
existence is well-known [19] . Numerical simulations of Luo-Hou [18] suggest very strongly 
that axisymmetric solutions can blow up: their model features initial data where both the 
vertical and angular velocities are odd in the vertical coordinate z, and they observe numer¬ 
ically a blowup at a fixed point on the boundary. The increasingly common view among 
experts [23111! is that functional analysis estimates are not sufficient to establish blowup, 
and instead one must analyze the geometry of trajectories in a careful way and use the spe¬ 
cial features of the Euler equation. We believe in addition that the Riemannian structure of 
the equation, as a geodesic equation on the group of volumorphisms as found by Arnold [lj, 
is also quite useful in both analyzing the local geometry and in finding simpler models with 
the same behavior. 

Our configuration manifold is the solid torus M = D x S * 1 = {(x, y, z) \ x 2 + y 2 < 1, —it < 
z < 7r}/{(x,y, —it) = (x,y,7i)}. The Euler equation for an axisymmetric velocity field U 
which is tangent to dM is given by 

(1) U t + U ■ VU = -VP, div U — 0, d e U = 0, (U, e r ) r=1 = 0, U t=0 = U 0 , 

where the pressure P is determined nonlocally by the equation 

(2) AP = — div (U ■ VU) = — Tr (VU) 2 , (VP, e r ) r=1 — —(U ■ VU, e r ) r=1 . 

The vorticity u = curl U satisfies the conservation law 

(3) + [U,u>\ = 0. 

It has been known since the work of Ebin and Marsden ra that the Euler equations have 
H s solutions for s > 5/2 on some time interval [0, T) for initial data Uq in H s . The most 
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famous criterion for extensibility of solutions past time T is the Beale-Kato-Majda result [3j 
that 

rT 


(4) 


\u 


LOO (It < OO. 


Expressing U in components as U = ae r + bee + ce z , we observe in Proposition |2.4| that 
if initially c is odd in z, a is even in z, and b is either even or odd in z, then the same 
remains true as long as the solution persists. In either case the points where z = 0 and 
either r = 0 or r = 1 become fixed points. Criteria for either blowup or global existence 
in a neighborhood of these fixed points have been studied by Chae 00 and by [7] in the 
case where b is odd, but to date we do not know of any other analysis of the case where b 
is even; as the latter case allows the vorticity to be nonzero at the fixed point, we believe it 
is quite interesting. For example, we prove the following theorem that severely restricts the 
possibilities for blowup at the origin when the angular velocity is even and nonzero. 


Theorem 1.1. Suppose Co is odd in z, and that ao and bo are both even in z. Then the 
origin is a fixed point of the flow (J5]) . Assume that the solution of ([Tj) ceases to exist at time 
t = T due to pointwise blowup of vorticity at the origin; that is, assume J^\u(t, 0)| dt = oo. 
Then the pressure function P given by (|2]) must satisfy limsup t _> T (T — t) 2 P rr (t, 0) > \. 


Superficially it seems easier to achieve blowup in the case where the angular velocity b is 
odd than when it is even, especially in light of Theorem 1.1 However the situation changes 
dramatically on the boundary r = 1 of the cylinder. 


In that case numerical evidence 
shows that the solution of (JTJ) ends in finite time due to blowup of odd data, yet our result 
below shows that blowup on the boundary is actually easier in the even case. 


Theorem 1.2. Let bi =6 0 (1,0) andbo = (6 0 ) r (l,0), and assume that 2bi(bi + b 2 ) = — c 2 < 0 
and a = a r (0,1, 0) < 0. Then either of the following is sufficient to ensure blowup: 

• —k 2 < P rr (t, 1, 0) + 3b 2 < 0, where k is the positive solution of k(k + a) = c 2 In 2, for 
all time; or 

• P r r(t, 1, 0) + 3b 2 > 0 for all time. 


It is also possible to prove blowup based on assumptions about the time derivative of 
pressure; although such assumptions are difficult to justify physically, it is interesting that 
relatively mild assumptions can have rather strong consequences. For example we show the 
following theorem. 


Theorem 1.3. Suppose Co is odd in z and ao is even in z, with b 0 either even or odd in z. 
Assume that P zz ( 0,1,0) < 0, that c z ( 0,1,0) < 0, and that c^(0,1, 0) 2 + P zz (0, 1, 0) > 0. If 
P zz (t , 1, 0) is increasing for all time, then the solution of ([TJ) ends in finite time. The same 
is true if ( t , 1 , 0) is replaced everywhere with (t, 0, 0). 

The motivation and context for the results described above comes primarily from the 
geometric point of view, in particular in terms of the local geometric behavior near a fixed 
particle trajectory. Generally, particle trajectories are given in terms of the Lagrangian flow 
p\ [0, T) x M —>• M by 

(5) Vt(t,x) = U(t,p(t,x)), rj(0,x)=x. 

The divergence-free condition div U = 0 implies that the Jacobian is preserved: 

(6) Jac(??) = det (Dp) = 1. 
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The vorticity transport equation ([3]) becomes, in Lagrangian form, the equation u (t, rj(t, x)) = 
Drj(t,x)uo(x), from which it is clear by Q that it is sufficient to understand when Dp ap¬ 
proaches infinity (or equivalently zero). Combining (|TJ) and (|5]) we obtain the Lagrangian 
form of the Euler equation, 


(7) 


D 


—= -VP(t,ri(t,x)), 


in terms of the covariant derivative. Differentiating ([7]) in space, we see that Dr / satisfies the 
linear ODE 


( 8 ) 


D D 


dtdt Dri ^ t,X ^ = Dr l(^ x )^ 


(The equation in this form is due to Ohkitani [21].) 

We now describe the Riemannian geometric interpretation of equation (J 8 ]) (based on the 
work of Arnold CD- The configuration space is the space of volume-preserving diffeomor- 
phisrns (volumorphisms) of M , a group under composition, with right-invariant L 2 Rie¬ 
mannian metric given by the kinetic energy; geodesics in this metric are given by curves rj 
satisfying the flow equation (|5]) with velocity field U solving the Euler equation ([Tj). The 
solution operator Uq r/(l) takes H s velocity fields to H s volumorphisms for s > §, and 
this map is called the Riemannian exponential map. Ebin-Marsden na showed that this 
map is C°°, and Ebin-Misiolek-Preston [13] showed that the map is not Fredholm, which 
means that its differential does not have finite-dimensional kernel or cokernel in general. 

The failure of Fredholmness is related to the possibility of infinite clusters of conjugate 
points, and in fact we see that accumulation of vorticity up to a blowup time T generically 
implies the existence of pairs t n ,t n+1 such that the geodesic 77 fails to minimize on [t n ,t n+ 1 ] 
for a sequence t n /* T. This should be viewed as a “positive curvature blowup,” since it 
cannot happen if the Riemannian sectional curvature is bounded above. In case this does 
not happen, we get fairly precise information on the growth of eigenvalues of the stretching 
matrix m■ 

Theorem 1.4. Suppose cq is odd and b 0 is even, with J^\u(t,r o ,0)\dt = 00 for either 
r 0 = 0 or r 0 = 1. Then there is an infinite increasing sequence t n /* T with p(t n ) conjugate 
to p(t n+ 1 ) if either 

• r 0 = 0 or 

• r 0 = 1 and 46 0 (1,0) 2 + 46 0 (1, 0)(6 0 ) r (l, 0) > |a rt (T, 1, 0)| 2 , 

where a is the radial component of the diffeomorphism p. In either case the supremum of 
sectional curvatures sup v K(U,V) approaches positive infinity in the sense that 

f T 

/ sup K(U(t), v) dt — 00 . 

Jo UeT id Diff M (M) 

Theorem 1.5. Suppose cq and h 0 are both odd, with Q |a>(t, r 0 ,0)| dt = 00 for either r 0 = 0 
or r 0 = 1. Then we cannot find curve-shortening variations on time intervals [a,b\ C [0,T] 
that are supported in arbitrarily small neighborhoods of the point r = r 0 , z = 0. In addition 
the pressure Laplacian and the sectional curvature must both approach negative infinity in 
the sense that 

(9) [ inf A P(t,x) dt — f inf K(U(t), V) dt — — 00 . 


UeT id Diff M (M) 
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We thus see that there is a basic dichotomy between positive-curvature blowup scenar¬ 
ios, corresponding intuitively to convex pressure and dominance of rotation, and negative- 
curvature scenarios, corresponding to concave pressure and dominance of stretching. 

We now describe the plan of the paper. First in Section [2j we present the axisymmetric 
3D Euler equation 0 in Lagrangian form. From here we compute the equation (| 8 ]) in 
components and specialize to the equations on the symmetry axis and on the boundary, 
where they simplify drastically. Then we discuss further simplifications that arise if the data 
is assumed to have additional reflection symmetries, in order to obtain a system of ODEs, and 
relate the resulting equations to well-known dynamical systems such as the Ermakov-Piimey 
equation. In Section [3] we collect conditions on the pressure Hessian leading to blowup at 
the fixed points and establish Theorems |I.1| 1 .3 [ Then in Section [4] we prove Theorems 


1.4 1.5, thus giving a geometric interpretation of our blowup results both in terms of the 


local geometry along trajectories and in terms of the Riemannian geometry of the volume¬ 
preserving diffeomorphism group. Lastly, in Section [5] we discuss future work by relating the 
axisymmetric 3D Euler equation to two lower-dimensional models we believe have similar 
geometric structure: In one dimension, an equation proposed by Wunsch as a geometric 
model of 3D Euler, based on an idea of Constantin-Lax-Majda [9] and developed by De 
Gregorio [10], and in two dimensions, the surface quasigeostrophic equation (SQG), which 
is used as a lower dimensional analogue of the 3D Euler equation. 


2. Background 


Let us first establish our notation. We denote the velocity held by U(t,r,z), and write it 
in components as U — ae r + beg + ce z , where e r , eg, and e z are the usual cylindrical unit 
vector fields and the coefficients a, b, c all depend on ( t,r,z). We denote the Lagrangian 
how r]: [0, T) x M -)• M by 

(10) rj(t,r,9,z)= (a(t,r,z),9 + fl(t,r,z),j(t,r,z)), 

with initial conditions a(0,r,z) = r, (3(0,r,z) = 0, and y(0,r, z) = z. The Lagrangian how 
equation ([5]) takes the form 

O' 

(11) cx t — a(t, a, 7 ), (3 t = -b(t,a, 7 ), 7 * = c(t, a, 7 ). 

r 

The following conservation law for angular vorticity is one of the most important for us. 


Proposition 2.1. For any axisymmetric solution U = ae r + beg + ce z of ([Tj) with flow 
components 0 , we have 

( 12 ) a(t, r, z)b(t, a(t, r, z), j(t, r, z)) = rb 0 (r , z). 

Proof. This comes from writing the angular component of (Jl| as b t + ab r + cb z + ab/r = 0. 
Composing with the how and using 0 . we get that dflaflt) = 0. Integrating and using 
a|t=o = r, we get □ 

Incorporating the conservation law of Proposition |2.1[ we get convenient forms for the a 
and 7 equations. 


Proposition 2.2. The components a and 7 of the flow (10) satisfy 

r 2 bl 


(13) 


Oft - r- — —P T (t. a, 7) and ~f tt = -P z (t, a, 7). 


a 
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while the Jacobian determinant condition 0 is 

(14) a(a r % - a z ^ r ) = r. 


Proof. The formulas (13) are proved by writing the equation (Jl| in components as a t + aa r + 
ca z — b 2 /r = —P r and c t + ac r + cc z = —P z , and using equations ( 11 ) to get second-order 
equations for g, then plugging in (12). The formula (14) is straightforward. □ 


There are effectively two boundaries on the cylinder: at r = 0 along the axis, and at 
r — 1. They are slightly different since on the axis the constraints are determined by 
the requirement of smoothness and rotational invariance, while on the boundary they are 
determined by the no-flow condition, although they are ultimately similar since a is fixed in 
either case: a(t, 0 , z) — 0 and a(t, 1 , z) = 1 . 


Lemma 2.3. At the fixed radii r = 0 and r 
the form 


1, the incompressibility condition (14) takes 


(15) a r (t, 0, z) 2 g z {t, 0, z) = 1 and a r (t, 1, z)g z (t, 1, z) = 1. 


Proof. This follows from the fact that a(t, r 0 , z) — e for r 0 = 0 or r 0 = 1, 
and (14) becomes (a/r)a r y 2 = 1. The term a/r becomes a r at r 0 = 
at r 0 = 1 . 


so that a z (t, e, z) = 0 
0 and becomes unity 
□ 


We obtain fixed points for the flow (for any 6) at z = 0 if we assume that the flow is 
symmetric about the plane z = 0 , in the sense that c 0 is odd in z\ that is, c 0 (r, —z) = 
—Co(r, z). The following result is well-known in the case when 5 0 is odd, but the case bo even 
has not been as widely studied. 


Proposition 2.4. Let U 0 be a divergence-free velocity field with components U 0 = a 0 e r + 
boee + coe z on the solid torus M = D x S 1 , and assume that ao is even, cq is odd, and b 0 
is either even or odd in the z variable (viewed as an element of [— 7 r, 7 r]j. Then the same 
symmetries hold for U for all time. In particular the origin r = 0 and z = 0 is a fixed point 
of the flow, as are all points on the circle z = 0 and r — 1 . 


Proof. Proposition 2.1 implies that if b q is even, then so is b 2 for all time. It is then easy 


to check that the equations (13) are preserved under the assumed symmetries, since writing 
out equation (| 2 ]) in components as 

(16) \d r (rP r ) + P zz = - a 2 /r 2 — a 2 - 2 a z c r — c 2 + 2 bb r /r 


shows that P will be even in 0 . We then have c(t, r, 0) = 0 for all t and r, so that by 0 
we have y(t, r, 0) = 0 for all t and r. Since we also have a(t, 0, z) — 0 and a(t, 1, z) = 1, we 
conclude that z — 0 and r = 0 or r = 1 are all fixed points. □ 


Note that the origin is a genuine fixed point, while the circle r = 1 and z = 0 is fixed by 
the flow but may rotate since b is not necessarily zero there (if it is not odd). 

Now differentiating the equations (13) spatially, we obtain linear ODEs for the components 


of the matrix Dp at these fixed points. 


Proposition 2.5. Suppose U 0 is a vector field satisfying the assumptions of Proposition \2fl 
Then on the axis f 0 (t ) = a r (t, 0, 0) and g 0 (t ) = 7 z (i, 0, 0) satisfy the ODEs 


(i7) m 


b 2 
11 o 


foflf 


= -P rr {t, 0, 0)/ 0 (t) and g'flt) = -P zz {t , 0, 0)^^) 
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with the constraints fo(t) 2 9o(t) = 1 and 

(18) 2P rr (t, 0, 0) + P zz (t, 0, 0) = —6 m 2 /f(t) 2 + 2bg//(t) 4 , 

where bo = <9 r 6o(0,0). And on the boundary, /i(t) = a r (t, 1,0) and g\{t) = 7 z (t, 1,0) satisfy 
the ODEs 

(19) f"(t) — 2bi(bi + b 2 ) + 3blf\(t) — —P rr (t, 1, 0)fi(t) and g'[(t) = -P zz (t, 1, 0)gi(t), 
with the constraints fi{t)gi(t) = 1 and 

P rr {t , 1,0) + P zz (t, 1, 0) = -3b? + 2bl( ^ b2) - 2 m 2 lh(t)\ 

where bi = 60(1,0) and b2 = (6 0 ) r (l,0). 


The equations (17) and (19) are essentially like the Jacobi equations in a three-dimensional 
manifold along a geodesic, where the Hessian of the pressure acts as the effective sectional 
curvature; this analogy is made more precise in [221 [23j. The equation © is essentially 
the Ermakov-Pinney equation [15] which describes a planar harmonic oscillator with central 
returning force. We have blowup if any of the terms fo(t) or go(t) become zero in Enite time 
(since the constraints mean that if one is zero, the other approaches infinity). Our main 
technique will thus be the comparison theory for linear ODEs. 


3. Blowup criteria 

Here we collect some conditions on the pressure Hessian which would ensure blowup at 
the fixed points of the equations. They are all based on the fact that some component of 
Dg must approach zero in order for the other component to approach infinity, and it is of 
course easier to ensure that the solution of a linear differential equation approaches zero than 
infinity since we need not require the coefficients to approach infinity. Geometrically we are 
asking that the Riemannian curvature is positive enough that we get conjugate points, since 
Di 7 generates the Jacobi fields and we want to see them vanishing. 

3.1. Blowup criteria on the axis. The following result generalizes the case e = 0 with 
P rr (t, 0,0) > 0 considered by Chae [5j; here we consider blowup in terms of a linear ODE 
for the component of Dr) rather than a Riccati equation for the velocity component. 

Theorem 3.1. Suppose c 0 is odd, a 0 is even, and b 0 is odd or even in z. Then any of the 
following is sufficient to ensure blowup in finite time, for either ro = 0 or ro — 1. 

• P rr (t, r 0 , 0) > 0 for all t, and initially a r (0,r 0 ,0) < 0 and 6 0 (0,r 0 ,0) = 0/ 

• P zz (t, r 0 , 0) > 0 for all t, and initially < 7 ,( 0 , r 0 , 0) < 0. 

Proof. All four cases are the same: if gift ) = a r (t, r 0 , 0) or g{t) = 7 z (t , r 0 ,0), the assumptions 
are equivalent to o(0) = 1, </(0) < 0, and q"(t)/q(t) < 0. Hence comparing with q 0 (t) = 
1 + g'(0)t we have g(T) = 0 no later than T = -l/ 5 '(0). □ 

Note that while Chae’s result [5] only requires P rr to be positive along the axis (not 
necessarily unbounded), equation ( [l 8 | ) clearly requires P zz to blow up to negative infinity in 
this case, corresponding to the fact that boundedness of the pressure Hessian is also sufficient 
to prevent blowup [Tj. 

In the case where (6 0 ) r (0,0) 7 ^ 0, the situation becomes more interesting. It is impossible 
in this case for a r to approach zero in finite time if P rr (t, 0,0) is bounded, and in fact 
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using fairly standard Sturm-Liouville comparison theory [26], we can get a blowup rate for 
P rr {t, 0,0) if a r does approach zero and satisfies a localized BKM-type criterion in the form 

f T f T dt 

/ |w(t, 0,0)| dt = 2\bo\ / = oo; 

Jo Jo 0, 0) 

see Theorem 14.21 _ 

Proof of Theorem |l.l[ Consider the equation p(t) — &o/p(£) 3 = — F(t)p(t ), where p(t) = 
a r (t,0,0) and F(t) = P rr (f,0,0). As pointed out by Eliezer and Gray [11], this is the 
equation for the radial coordinate p(t ) for a planar central force system x(t ) = —F(t)x(t), 
y(t ) = —F(t)y(t), where pit) 2 = x{t ) 2 + y(t ) 2 and x(0)y(0) — x(0)r/(0) = b 0 . The angular 
coordinate 9{t ) then satisfies 9{t ) = b 0 /p(t ) 2 , and the condition J Q T dt/p{t) 2 = oo means 
exactly that 9{t ) winds around the origin infinitely many times as t T. Hence x(t) and 
y(t) have infinitely many zeroes as t —>• T, and thus every solution of g(t ) = — F(t)g(t ) has 
infinitely many zeroes on (0, T). 

We now change variables so that the blowup time is sent to infinity: set s = — In (T — t), 
set g(t) = (T — t) 1//2 j(s), and F(t) = H(s)/(T — t) 2 . Then the equation g{t) = — F{t)g(t ) 
becomes 

i"0) = Q ~H{s)) j(s), 

and we must have 1 iin sup,^^ H(s) — ^ > 0 for solutions of this equation to have infinitely 
many zeroes. Translating back in terms of P rr we obtain 

(20) limsup(T — t) 2 P rr (t : 0) > |. 

t^T 


□ 

With more assumptions one can obtain more precise criteria, using the methods presented 
e.g., in Swanson [26], but for our purposes (20) already makes clear how tightly the pressure 


Hessian is constrained in a typical blowup scenario at a fixed point with nonzero vorticity. 

Theorem L2 shows that we can get the same sort of blowup as in Theorem 3T even if the 
pressure Hessian is negative, as long as there is enough rotation to compensate. 

First we work under the condition — k 2 < P rr (t , 1,0) + 3b 2 < 0. 


1.2 


Proof of Theorem_ 

Write fi{t) = a r (t, 1, 0) and F(t) = P rr (t, 1, 0) + 3bf; then by equation (19), fi satisfies the 
equation 

(21) ff{t) = -c 2 - F(t)flit), hi 0) = 1, /](0) = -a, 

where 0 < —Fit) < k 2 and a > 0. Consider the solution yi(t) of the related problem 

Vi(t) = -F(t)yiit), (0) = 1, j/i(0) = 0. 

By the usual Sturm comparison theorem [26], we have 


( 22 ) 


1 < i/i (t) < cosh kt for all t. 


Using the reduction of order trick, the solution of 1/2 (t) = —Fit)y 2 {t) with t/ 2 (0) = 0 and 
2/2(0) = 1 is given by 2/2(^) = Vi(t) Jq ds/yiis) 2 . Using variation of parameters, we may write 
the solution h(t) of (21) as 

flit) = Vi(t) ~ ay 2 (t) + c 2 yiit) [ y 2 (s) ds - c 2 y 2 it) [ yi(s)ds 


= Vi(t) 


1 — a 


ds 


t rt 


yi(s) 


— c 



yijs) 

yiir ) 2 


dr ds 
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The inequality (|22|) now implies that 

f* ds 


h(t) 

2/i (*) 


< 1 


'0 


cosh' ks 



dr 


= 1 


a , , 

— tanh kt 

k 


■ :1 In 2 c 


0 Js 
,2 


cosh' kr 


ds 


k 2 


+ F ln(l + e - 2 “) + 


2 cH 


k( 1 + e 2kt ) 


and thus we have 


fi (t) k 2 — ak — c 2 In 2 
Inn —t < - - 2 -. 


t^oo yi (t) 

We conclude that fi(t) is eventually negative. 

The second assumption that P rr (t, 1, 0) + 3b 2 > 0 is much easier: in this case we just have 
1 , 0 ) < —c 2 , and since ct r ( 0 , 1 , 0 ) = 1 and a tr ( 0 , 1 , 0 ) < 0 , we obviously have a r (t, 1 , 0 ) 
reaching zero in finite time. □ 

In the previous theorems, we showed how some assumptions on the sign of certain com¬ 
ponents of the pressure Hessian could cause blowup on either the axis or on the boundary. 
We now change our perspective somewhat and look at the consequences of assumptions on 
the time derivative of components of the pressure Hessian. 

Proof of Theorem |1.3| . From Proposition |2.5| we see that the simplest case is the y 2 
component, since in all cases (b even or odd, or whether we work on the boundary or the 
axis), 7 z satisfies the simplest equation y tt2 = — P zz (t)lz- Hence we need only consider the 
equation 

(23) g"{t) = -Q(t)g(t), g{ 0 ) = 1 , g'{ 0 ) = -a, 

where a > 0, Q(0) < 0, and v 2 = a 2 + Q(0) > 0, and we assume Q\t ) > 0 for all t. 
Multiplying (23) by g' and integrating, we obtain 


g\t? + Q(t)g(t) 2 = h 2 + / Q'{r)g{r) 2 dr. 

Jo 

Rearranging this now gives 

g'(t) 2 -V 2 = -Q{0)g(t) 2 + f Q'{t) [g{rf - g(t ) 2 ] dr. 

Jo 

Since g(0) = 1 and (/(0) < 0, we know g is decreasing and positive on some time interval 
[0, T 0 ]. On this interval we have g{r) 2 — g(t ) 2 > 0 whenever r < t < T 0 . We conclude that 
on this interval, g'{t) 2 — v 2 > 0, which implies that g'{t ) < —v as long as g(t) is decreasing 
and positive. Since g' is continuous, we conclude that g must reach zero before it changes 
direction, and furthermore we have g(t ) < 1 — vt so that the time T of the first zero is no 
larger than 1/v. □ 

Of course the same theorem applies with a r replacing y 2 , in case bo is assumed to be 
odd, since then the equations all reduce to the same equation f"(t ) = — Q(t)f(t ). 

Similar theorems could be proved using the more complicated equations for a r arising from 


(17) and (19) in case b 0 is even, but we will leave these aside for now. 


4. Global geometry of the Euler equation 


In the previous Section we made a variety of assumptions on the local behavior of the 
fluid which could lead to blowup; here we would like to tie this local picture into the global 
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behavior of the equation (especially as related to the Riemannian geometry of the volume¬ 
preserving diffeomorphism group) and the global behavior of the pressure function. 


4.1. Conjugate points and blowup in axisymmetric fluids. Viewed as a Riemann¬ 
ian manifold, the group of volumorphisms Diff M (M 3 ) with Riemannian metric {{u,u)) = 
Im (u,u) dp, has geodesics satisfying (|7j) (which is equivalent to (J5]) and ([TJ)) , as pointed out 
by Arnold [Tj. Its sectional curvature describes small Lagrangian perturbations, and is given 
for divergence-free velocity fields U and V, using the Gauss-Codazzi formula j2Q], by 


(24) 


((R(U,V)V,U)) L 2= / V 2 P(V,V)dv~ / |VQ| 2 dp, 


IM 


' M 


where P is the pressure of the velocity held U and Q solves the Neumann problem A Q = 
— div (U ■ W) with (VQ + U ■ W, e r ) r =i = 0. In the sense that the Hessian of the pressure 
always plays the role of the “effective” curvature along local Jacobi fields, as mentioned after 
Proposition |2.5[ the formula (24) shows that the actual Riemannian curvature is always less 
than this “effective” curvature term. 

The following theorem shows that conjugate points on the volumorphism group can be 
found by a local criterion along a single particle path, which demonstrates that it is easy 
to shorten a geodesic curve by performing rotational perturbations near a point. This is a 
purely three-dimensional result: Fredholmness of the 2D Riemannian exponential map m 
implies there is no similar result in two dimensions. 


Theorem 4.1. [22j Suppose M is a three-dimensional manifold and x is in the interior of 
M, and p is a Lagrangian solution of the ideal Euler equations with rj t (0,x) = uq{x). Let 
A [t,x) = Drj(t,x)lDgft, x), and let uio(x) = curlt/ 0 (a;) denote the initial vorticity at x. If 
there is a vector field vft) along the Lagrangian trajectory 1 1 -> p(t,x) with v(fi) = v(t 2 ) = 0 
such that 

rt2 

(25) I(v,v ) = / (A (t, x)v(t),v(t)) + (cuo(a;) x v(t),v(t)} dt < 0, 

Jt i 

then rj(ti) is conjugate to rj(r) for some r <t 2 ; in particular the geodesic rj is not minimizing 
on [ti,t 2 ]. 


In general for axisymmetric flows in the orthonormal basis {e r , eg, e z } we have 

/ a r 0 a z \ 

Drj(t,x ) = a/3 r a/r a/3 z , 

\ 7r 0 ) 

and since ^ for f3 reduces to |j^( t,r,z ) = b 0 (r,z) by the conservation law (12), we have 
j3 r (t,r,z) = t(h 0 ) r (r,z ) and (3 z (t,r,z ) = t(b 0 ) z (r, z). Thus we have 


a r 0 0 


a r 0 0 

Z>//(f,l,0) = ( t(6 0 )r(l,0) 1 t(6 0 ) 2 (l,0) ) . 

0 0 7 , 


(26) £>77(t,0,0) = ( 0 a r 0 | , 

0 0 

We now analyze the index form on the axis and on the boundary. In all cases we assume 
the localized Beale-Kato-Majda criterion (31): at the fixed point z — 0 with either r 0 = 0 or 
ro = 1, we have 


| u>(t, Tj(t, r 0 , 0)) | dt = oo for r 0 = 0 or r 0 = 1. 
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Note that if bo is odd, then the vorticity is identically zero at the origin, and the localized 
BKM criterion cannot be satisfied there. We first consider the assumptions necessary to get 
an infinite sequence of conjugate pairs (and thus sectional curvature increasing to positive 
infinity). 

Proof of Theorem 11.41 First consider the situation at r = 0. It is sufficient to show that 
for any t\ > 0 there is a t 2 > t\ such that the index form I(v, v) in (25) can be made negative 
for some v vanishing at both t\ and t 2 - The initial vorticity is given by ujq = 2boe z , and 
since it is stretched by u(t, 0 , 0 ) = 7 z (t, 0 , 0 )u;o( 0 , 0 ) = 2bo/a r (t, 0 , 0 ) 2 , our assumption yields 
Jq dt/a r (t, 0, 0 ) 2 = 00 . Equation (26) yields 

A(£, 0,0) = Drj(t, 0, 0)' Drj(t, 0, 0) = ( 0 0 % 0 ) , 

V 0 0 7z V 


so that the index form (25) becomes, for v(t) = f(t)e r + g(t)eg + h(t)e z , 
I(v,v)= I a 2 J 2 + a 2 r g 2 + ^ 2 h 2 + b 0 (fg - gf)dt. 


H 1 


Set h = 0, integrate by parts using v(ti) = vfo) = 0, and complete the square to obtain 


f*t2 


I(v,v) = I (a r g+^f) z + a 2 f 2 -^f 2 dt 


We choose g = k — bo//a 2 , where k is chosen so that f t 2 g dt = 0, and obtain 


(27) 


I(v,v) = 


b 2 
11 0 


f*t2 


rt 2 


(*2 - «l) 2 Jt, 




For the latter integral, we rescale our time variable by s(f) — / 0 ( dr/ay(r) 2 , and obtain 


(28) 




a rf 2 ~ a^f 2 dt = I f\s) 2 -b 2 0 f( S ) 2 ds. 


r>S 2 


' 11 


By assumption we have s(t) —> 00 as t —y T, so that we can certainly choose t large enough 
so that with f(s) = sin2(.s — s 1 ) 7 r/bo, the integral in (28) is negative. With this choice, the 


a 2 + b 2 t 2 b 2 t 0 
b 2 t 1 0 

0 0 7 2 


Erst integral in (27) vanishes, and we get I(v,v ) < 0. 

Next we consider the situation at r = 1. We then have 

A(M,o) = 

and the initial vorticity is o; 0 (l,0) = (bi + b 2 )e z . For v(t) = f(t)e r + g(t)eg + h(t)e z , the 
index form (25) becomes 

I(v,v)= [ (a 2 + b\t 2 )f 2 + 2b 2 tf g + g 2 + 7 ?/i 2 + (bi + b 2 )(fg — gf) dt. 

Jt\ 

Set h — 0, integrate by parts, and complete the square to obtain 


(29) 


r*2 


I{v,v)= (g + (bi + b 2 )f + b 2 tf) + a 2 / 2 - bi(bi + b 2 )/ 2 dt. 


Hi 
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We have u>(t, 1,0) = e z so that the blowup condition is f^dt/a r (t, 1,0) = oo. If 

q := a rt (T, 1, 0) ^ 0 then the dominant term in (29) looks, for t\ and t 2 sufficiently close to 
T, like 

/' ! t 2 (T-trm 2 -im 2 dt 


for ( 2 = 4bi(bi + b 2 ). Minimizers of this integral subject to f(t\) = /(f 2 ) = 0 satisfy the 
equation A(g 2 (T — t) 2 f(t )) + ^/(t) = 0 , with solutions 


(30) 


/(*) = 


VT^t 


cos In (T — t) + 0 ) for 0 = \JC, 2 jq 2 — 1 


and some constant 0, and all such solutions vanish infinitely many times up to time T. 

As in the case above when r = 0, we choose g = k — b\f — b 2 4 (tf) where k is chosen 


so that j t] 2 g{t) dt = 0 , i.e., k = ^ f(t)dt. 


Then the first term in (29) vanishes if 


f t 2 f{t ) dt = 0, and we can easily choose functions / of the form (30) with vanishing mean 


if t 2 is close enough to T. 

The result on the integral curvature follows from rather general principles of compar¬ 
ison theory in ODEs; the following is adapted from the proof of Theorem 5.1 in Chap¬ 
ter 11 of Hartman [16] . Consider a Jacobi field J(t) satisfying the equation J"(t) = 
—i?( 7 (t), J(t)) 7 (f). Let 0(t) = || J(t)| 2 ; then we have 


0"(t) = - 2 K{t)<f>(t) + \J'(t)\ 2 > -2K(t)<f>(t), 


assuming that 7 (t) is a unit vector, where K(t) is the sectional curvature in directions J(t) 
and 7 (t). On an interval [t n ,t n+ 1 ] where a Jacobi field vanishes, we have 

ft ftn +1 

(W 1 - t n )(j>(t) <2 (t n+ 1 - t)(s - t n )K(s)<j>(s ) ds + 2 (t - t n )(t n+ 1 - s)K(s)<j>(s) ds 

•J tfi J t 

for any t G [t n ,t n+ 1 ]. If t 0 denotes the location of the maximum of 0(t) then using t — t 0 
and overestimating the right side by 0 (s) < 0 (t o ), then cancelling 0 (t o ), we get 

ftn+l 

(t n+ 1 - t n ) < \{t n+ 1 - t „) 2 / AT(s) ds. 

»/ t n 

Therefore we have 


^n + 1 


K (s) ds > / K (s) ds > 


n=l 


n=l 


'n+l 


Since the sum — tn) < T converges, the reciprocal sum must diverge, and we 

conclude that f 0 K(t) dt is positive infinity. □ 

More generally we can try to generalize the preceding computations to any Lagrangian 
path along which the vorticity is increasing to infinity, but it turns out there are not always 
infinitely many conjugate pairs along the trajectory. The alternate possibility is however 
rare enough that we can get fairly concrete information about the growth of the components 
of the stretching matrix, as in the following theorem of the first author. 
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Theorem 4.2. [24J Assume that a solution of the 3D Euler equation |l| on a 3D manifold 
M has a maximal existence time T < oo, and that the following strong form of the Beale- 
Kato-Majda criterion Q holds: 

(31) zte G M\dM s.t. f \ojft, rjft, x)) | dt — oo. 

Jo 

Then either there is a sequence t n of times such that rj(t n ) is conjugate to r]ft n+ i ) for every 
n (in other words, p fails to be a locally minimizing geodesic on the interval [t n ,t n+ i]), 
or there is a basis {ei,e2,e3} such that ujo(x) is parallel to 63 and the components A^ of 
Aft) = Dpft,x)^Drift, x) satisfy 


(32) 


A33 ft) 


/o A n W+A 22 (t) 


dt < 00 


and 


lim IlTDlA = 0 . 

t ~’ T Jo A 22 (r)dr 


Conjugate points imply that the curvature is approaching positive infinity, while the alter¬ 
native condition allows for negative curvature but implies that the stretching matrix A must 
have its eigenvectors aligning in fixed directions rather than rapidly rotating. In the present 
context we can get more explicit information along our trajectories due to the symmetry. 
Proof of Theorem |1.5[ The cases where r 0 = 0 and r 0 = 1 are similar, so we will just deal 
with r 0 = 1 here. If bo is odd, then we have by (26) that 

.,2 


A(t,i, OH 


0 0 

1 tb 3 

f2t-2 


tb 3 t 2 b'i + 7 2 

with the initial vorticity cuo(l,0) = b3e r , where 63 = (&o)z(l,0). Plugging vft) = f(t)e r + 
gft)ee + h(t)e z into the index form (25), we get 

r £2 

I(v,v) — / ay/ 2 + (g + tb 3 h + b 3 h)~ + y 2 /} 2 dt, 

Jt 1 

which is always positive. Hence if the Jacobi held is spatially supported in a sufficiently small 
neighborhood of the circle r = 1 and z — 0 , its index form in the volumorphism group will 
also be positive, and thus we cannot have conjugate points arising from locally-supported 
Jacobi fields. 

By formula (16), using the fact that P r (t, 1, 0) = bft, 1, 0) 2 = 0 by the boundary condition 
and aft, 1,0) = a z ft, 1,0) = bft, 1,0) = 0 by our oddness assumptions and no-how at the 
boundary, we get 

a tr (t, 1, 0) 2 


(33) 


A P{t, 1,0) = P rr ft, 1,0) + P zz ft, 1, 0) = -2 a r ft, 1, 0) 2 = -2 


a r (t, 1, 0) 2 

Since the vorticity at r = 1 and z = 0 is ojft, 1,0) = — b z ft, 1, 0)e r = — b 3 a r (t, 1, 0)e r , our 
assumption on the localized Beale-Kato-Majda criterion becomes J Q T a r (t, 1,0) dt = 00. Thus 
by equation (33) we have 

[ AP(t,l,0)dt = -2\\f(t)\\ 2 L 2 [0T] < -C sup |/(t)| 2 , 

Jo 0 <t<T 


where /(t) = In a r (t, 1, 0), with the last inequality following from the Poincare inequality on 
[0, T], Hence we must have f 0 A Pft, 1, 0) dt = —00 or else f(t) would be bounded, which is 
impossible since J Q T e^O dt = cx) by assumption. 
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To show that the curvature also approaches negative infinity, it is enough to show that we 
can find a divergence-free axisymmetric velocity field V such that K(U, V) < CAP(t, 1, 0) for 
some C independent of time. To do this we note that K(U, V) <C f M V 2 P(V, V) dfj,/\\V\\ 2 L2 , M j 
by formula (24), for some C which depends on the L 2 norm of U (which is constant in time). 
Thus it is enough to make f M V 2 P(V, V ) dp close to AP(1, 0). Let ( be a smooth function 
on M supported in (—1,1) which is odd, and for a small £ > 0, set if)(r,z) = C(~)C(~)- 


is ax- 


Then ijj is odd in z and vanishes when r = 1, so that the vector field V = —^fe r + 
isymmetric, divergence-free, with vertical component odd in z and tangent to the boundary 
r — 1. It is easy to check that 


'M 


V 2 P(V, V)dn=n c \nf duj ( f ((v) 2 dv'j (p rr (l, 0) + P zz ( 1, 0)) + O(e), 


which is what we needed to show. Hence formula (J9]) follows. □ 

Roughly speaking, if the swirl component b is even, then localized blowup may be associ¬ 
ated to extreme positive sectional curvature on the diffeomorphism group, while if b is odd, 
then it must be associated with negative curvature. Note that the condition bi(bi + b 2 ) > 0 
is the same condition b(r)oj(r) > 0 at r = 1 as the condition in [22J for an axisymmetric fluid 
flow u = b(r)eg to have positive sectional curvature in all directions, so that it is natural 
that this condition yields infinitely many conjugate points. Curiously this condition is also 


the opposite of the condition bi(bi + b 2 ) < 0 in Theorem 1.2 which ensures blowup 


5. Outlook 


We have seen that there are several local criteria for the pressure Hessian that would 
ensure a “local blowup,” which is to say a blowup of the vorticity along a single particle 
trajectory. Further we have analyzed the connection between this possible blowup and 
the index form along the geodesic in the volumorphism group, inspired by the fact that 
the index form involves precisely the quantities of vorticity and stretching which appear in 
the Beale-Kato-Majda criterion for blowup. The reason the index form along the actual 
geodesic in the volumorphism group can be approximated by a local index form along a 
single particle trajectory, and thus related to the supremum of vorticity, is precisely the lack 
of Fredholmness. 

As such we would like to analyze other PDEs which happen to represent geodesic equations 
on infinite-dimensional groups, which share the properties that their Riemannian exponential 
maps are smooth but not Fredholm. Such PDEs would be expected to be good models of 
the phenomena we observe here for the axisymmetric 3D Euler equation. To date there 
are only two other PDEs with these properties. In one dimension we have the Wunsch 
equation for which all solutions blow up in finite time iza due to what appears to 

be positive sectional curvature everywhere. In two dimensions the surface quasi-geostrophic 
equation (SQG) has a smooth exponential map which is not Fredholm, a recent result of 
Washabaugh [2S]. No global existence or blowup result is known for this equation, and even 
the analogue of Theorem |4.2 is not yet known, though it seems likely to be true. In fact all 
three equations are likely to have very similar geometric structure. 
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